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Abstract 

We present a subtraction scheme for computing jet cross sections in electron- 
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1 Introduction 

In recent years a lot of effort has been devoted to the extension of the subtraction method 
of computing QCD corrections at the next-to-leading order (NLO) accuracy to the com- 
putation of the radiative corrections at the next-to-next-to-leading order (NNLO) [H El El 
IH 13 El HI El E| • In particular, in Ref. [TU], a subtraction scheme was defined for computing 
NNLO corrections to QCD jet cross sections to processes without coloured partons in the 
initial state and arbitrary number of massless particles (coloured or colourless) in the final 
state. That scheme can be summarized as follows. 

The NNLO correction to any m-jet cross section is a sum of three contributions, the 
doubly-real, the one-loop singly-unresolved real-virtual and the two-loop doubly-virtual 
terms, 



Here the notation for the integrals indicate that the double-real corrections involve the fully- 
differential cross section d<7^+ 2 of m + 2 final-state partons, the real-virtual contribution 
involves the fully-differential cross section for the production of m + 1 final-state partons at 
one-loop and the double virtual term is an integral of the fully-differential cross section for 
the production of m final-state partons at two-loops over the phase space of m partons. The 
phase spaces are restricted by the corresponding jet functions J n that define the physical 
quantity. 

In d — 4 dimensions the three contributions in Eq. (jl.ljl are separately divergent, but 
their sum is finite for infrared-safe observables. (The requirement of infrared safety implies 
certain analytic properties of the jet functions J n that are spelled out in Ref. jH].) As 
explained in Ref. fU] we first continue analitically all integrals to d = 4 — 2e deminesions 
and then rewrite Eq. (jl.lj) as 




(1.1) 




(1.2) 
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that is a sum of integrals, 
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each integrable in four dimesions by construction. Here da^'^ 1 and da^'^ 2 are approxi- 
mate cross sections that regularize the doubly-real emission cross section in the one- and 
two-parton infrared regions of the phase space, respectively. The double subtraction due 
to the overlap of these two terms is compensated by da^'^ 12 . These terms are defined 
in Ref. ^U] explicitly, where the finiteness of da^fif is demonstrated also numerically for 
the case of e + e~ — *3 jets (m = 3). In Ref. [TT], we computed the integral J x da^'^ 1 and 
showed that the terms in the first bracket in Eq. (jl.4|) do not contain e poles. Neverthe- 
less, those terms still lead to divergent integrals due to kinematical singularities in the 
one-parton unresolved parts of the phase space. In this paper we define explicitly da^^ 1 

and ( f ± da^^ 1 ^) A \ that regularize the singly-unresolved limits of the real-vitrual cross 
section and j x da^'^ 1 in turn. Thus we complete the presention of all formulae relevant 

.NNLO 
m+l 



for constructing da^+i° explicitly. 



2 Notation 



2.1 Matrix elements 



We consider processes with coloured particles (partons) in the final states, while the initial- 
state particles are colourless (tipically electron-positron annihilation into hadrons). Any 
number of additional non-coloured final-state particles is allowed, too, but they will be 
suppressed in the notation. Resolved partons in the final state are labelled by i,k, I, ... , 
the unresolved one is denoted by r. 

We adopt the colour- and spin-state notation of Ref. [T2]. In this notation the am- 
plitude for a scattering process involving the final-state momenta {p}, \M m ({p})) , is an 
abstract vector in colour and spin space, and its normalization is fixed such that the squared 
amplitude summed over colours and spins is 

\M m \ 2 = (M m \\MJ . (2.1) 
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This matrix element has the following formal loop expansion: 

\M) = \M (fS >) + \M ( > 1) ) + ... , (2.2) 

where \A4^) denotes the tree-level contribution, is the one-loop contribution and 

the dots stand for higher-loop contributions, which are not used in this paper. 

Colour interactions at the QCD vertices are represented by associating colour charges 
Tj with the emission of a gluon from each parton i. In the colour-state notation, each 
vector \M.) is a colour- singlet state, so colour conservation is simply 

(T,T^\M) = O t (2.3) 

^3 

where the sum over j extends over all the external partons of the state vector \M), and 
the equation is valid order by order in the loop expansion of Eq. (|2.2|) . 

Using the colour-state notation, we define the two-parton colour-correlated squared tree 
amplitudes as 

\M^ k) {{p})\ 2 = (M (0) ({p})l TvT k \M®({p})) (2.4) 

and similarly the three-parton colour-correlated squared tree amplitudes, l-M^n] 2 for i, k 
and I being different, and the doubly two-parton colour-correlated squared tree amplitudes 
\M {0) I 2 - 

,(o) 



and 



\M$ km \ 2 = (-M (0) | {T 4 ■ T k ,T, ■ Ti} \M®) , (2.6) 



where the anticommutator {Ti-T k , Tj-Ti} is non-trivial only if i — j or k — I, see Eq. (J2.8I) . 
We shall also use the two-parton colour-correlated one-loop amplitude, defined using an 
analogous notation: 

2Re(M [0] \\M {1) ) m = 2Re(M i0) \Ti-T k \M {1) ) . (2.7) 



The colour-charge algebra for the product (Ti) n (T k ) n = Ti-T k is: 

T l -T k = T k -T l if i^k- T? = a- (2.8) 

Here C\ is the quadratic Casimir operator in the representation of particle % and we have 
C F = T R (N? - 1)/N C = (iV c 2 - 1)/(2N C ) in the fundamental and C A = 2T R N C = N c in the 
adjoint representation, i.e. we are using the customary normalization Tr = 1/2. 
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2.2 Dimensional regularization, one-loop amplitudes and renor- 
malization 



We employ conventional dimensional regularization (CDR) in d = A — 2e space-time dimen- 
sions to regulate both the IR and UV divergences, when quarks (spin-| Dirac fermions) 
possess 2 spin polarizations, gluons have d — 2 helicity states and all particle momenta are 
taken as d- dimensional. 

Turning to the renormalization of the amplitudes, let the perturbative expansion of the 
scattering amplitude A m in terms of the bare coupling g s = y/A~KQ^ be 



An 



I "4-771 ) 



</i 2£ 

An 



(2.9) 



where q is a non-negative integer and fi is the dimensional-regularization scale. For the 
renormalized amplitudes (in the CDR scheme) we use the notation \Ai m ). These are 
obtained from the unrenormalized amplitudes by expressing the bare coupling in terms of 
the running coupling a s (yU.|j) evaluated at the arbitrary renormalization scale as 



«s(/4) 

Air 



y + 0K 2 (/4)) 



(2.10) 



where (3q is the first coefficient of the j3 function for rif number of light quark flavours, 

(2.11) 



Po = y C A - ^T R n { - ^T R n s , 



for QCD n s = 0. In Eq. ()2.10|) . S E is the phase space factor due to the integral over the 
(d — 3) -dimensional solid angle, which is included in the definition of the running coupling 
in the MS renormalization scheme* 



S d -Vfl (4tt) 



d-3 



ra-e) 



(2.12) 



We always consider the running coupling in the MS scheme defined with the inclusion of 
this phase space factor. 

The relations between the renormalized amplitudes of Eq. (J2.2j) and the unrenormalized 



*The MS renormalization scheme as often employed in the literature uses S e = (47r) e e _e7E . It is not 
difficult to check that the two definitions lead to the same expressions in a computation at the NLO 
accuracy. At NNLO these lead to slightly different bookkeeping of the IR and UV poles at intermediate 
steps of the computation, but the physical cross section of infrared-safe observables is the same. Our 
definition leads to somewhat simpler bookkeeping at the NNLO level. 
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ones are given as follows: 



' ( 2 \ 2e \ <?/ 2 

|M<?> = ( %^S e -M |^>, (2.13) 



Itt 



\M$) = {^J^S^Y (ft \A%) - q -^S £ |^>) • (2-14) 

After UV renormalization, the dependence on \x turns into a dependence on fi R , so 
physical cross sections depend only on the renormalization scale fi R . To avoid a cumbersome 
notation, we therefore set fi R = \x in the rest of the paper. Furthermore, after the IR poles 
are canceled in an NLO, or NNLO computation, we may set e — 0, therefore, the fi 2 R and 
S' 1 factors that accompany the running coupling in the renormalized amplitude do not 
give any contribution, so we may perform the 

1% s l Y" (".(tids-iV _ f^M] q/2+l (2.15) 



47T 6 J V 4tt ~ f - J \ 4ir 
substitution in Eqs. (l2~T3j) and (|2~T3J) . 



2.3 Remark on regularization-scheme dependence 

Although the application of conventional dimensional regularization (CDR) is conceptually 
clean, the computation of squared matrix elements is much simpler in other versions of 
dimensional regularization, most notably in dimensional reduction (DR). As a result, most 
of the multiparton QCD amplitudes |-4m' ) ), both at tree-level (n = 0) and one-loop (n = 
1), are available in DR. At the level of cross sections however, the CDR scheme is used 
traditionally, therefore, the relation between the two schemes has to be established. 

The regularization-scheme (RS) dependence of the matrix elements at tree level affect 
only terms of 0(e), therefore, in computing the (m + 2)-parton cross section in Eq. (jl.3j) 
the RS dependence is completely harmless, the difference vanishes when we take the four- 
dimensional limit. The subtraction terms that regularize the real emission also depend 
on the RS. While this dependence does not influence d<7^+2°> ^ leads to differences (even 
in divergent terms) when the subtraction terms are integrated over the factorized phase 
space of the unresolved parton(s). The standard practice in the literature is to set up the 
subtraction scheme in CDR and transform the loop matrix elements to CDR if those were 
obtained in other schemes, for instance, in DR. 

The RS dependence in the loop amplitudes has in general both ultraviolet (UV) and 
infrared (IR) origin. Both have been discussed thoroughly up to two loops in Ref. jTHj. In 
the present paper we deal only with one-loop amplitudes and we summarize the transition 
rules from DR to CDR here. 
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The UV part of the RS dependence is due to the RS dependence of the renormalization 
procedure. At the one-loop level it means that Eq. (J2.14)) remains valid, with the same 
expansion parameter, no matter in which RS the bare amplitudes are computed if we 
perform the substitution [T2| 

Po^Po + effi 3 (2.16) 

in Eqs. (J2~TUJ) and (f2~Ti|) . By definition in CDR /3^ DR = 0. If the bare amplitudes are 
computed in the DR scheme, then /3° R = — Ca/6. 

The IR part of the RS dependence can be decomposed into universal finite terms and 
non- universal contributions at 0(e) ^3]. The finite terms are completely factorized, 

I^S£(a* 2 ; M)> = \ (S>* RS ) I<rs(^ 2 ; M)> + l*£V; M)> + o(e) , (2.17) 

while the 0(e) contributions do not contribute to dcr^ 1 ^ in the four-dimensional limit. 
The transition coefficients that relate the amplitudes in the RS's depend only on the flavour 
of the external partons and were first computed in Ref. [T^]. If 7j CDR = 0, as always assumed 
by definition, then 

~DR _ ~DR _ 
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2.4 Cross sections 



% DR = ^, 7 9 DR =^- (2-18) 



In our notation the real-virtual cross section dcr R ^ 1 is given by 

daZi = #m+i(M) 2Re(M < £ +1 \M < £ +1 ) , (2.19) 

where d(j) m +i({p}) is the ^-dimensional phase space for m + 1 outgoing particles with 
momenta {p} = {pi, . . . ,p m +i} an d total momentum Q, 



Q) 



m+1 d d p- 

n (2#t 5 +^ 



(2TT) d 5^( Pl + ---+p m+1 -Q). (2.20) 



The integral of the singly-unresolved approximate cross section for doubly-real emission 
over the factorized one-parton phase space was computed in Ref. 

^da^ 2 Al =da R +1 ®/(m,e), (2.21) 

where d<7 R +1 is the Born-level cross section for the emission of m + 1 partons and I(m, e) 
is an operator acting on the colour space of the m + 1 final-state partons. The notation on 
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the right hand side means that one has to write down the expression for dcx^ +1 and then 
replace the Born-level squared matrix element 



\M 



(0) |2 
m+ll 



(0) 



+1/ ' 



(2.22) 



by 



(M { Zi\I(m,e)\MZi)- 



(0) 



The insertion operator I(m, s) differs from the 1(e) operator derived in Ref. 
singular terms as e tends to zero. Explicitly, 



(2.23) 
in non- 



I({p};m,e) = ^S e [%2 



2 \ £ 



X 



Ci(y iQ ; m, e) T\ + ^ S ik (y ik , y iQ , y kQ \ m, e) T{T k 



(2.24) 



where y ik = s ik /Q 2 = 2pi-p k /Q 2 , y iQ = 2pi-Q/Q 2 and with 



— C qg — CS 



C g — ^Cgg + n (Cqq ~ CS . 



(2.25) 



Explicit expressions for the functions C ik (y iQ ;m,e), S ik (y ik , y iQ , y kQ ; m, e) and CS(m,e), 
can be found in Ref. where it was shown that the e poles of the one-loop squared 



matrix element 2 Re(A^^l l _ 1 |A^^+ 1 ) are cancelled exactly by (M^' +1 \I({p};m,e)\A4^ +1 ) 



(0) 



,(0) 



3 Counterterms for the real-virtual cross section 
3.1 Factorization in the collinear and soft limits 

In order to devise the approximate cross section dcr^+'f' 1 , we have to study the factorization 
properties of one-loop squared matrix elements when one parton becomes soft or collinear 
to another parton. The relevant factorization formulae have been computed in Refs. |161 
El EU El 121] • In our work we use the formulae of Ref. [20] for collinear parton splitting 
and those in Ref. j2I] for soft gluon emission. However, the notation in those papers is 
not convenient for writing factorization formulae which avoid double counting in the soft- 
collinear limit, therefore, we present new formulae here. 



S 



3.1.1 Collinear limit 



We define the collinear limit of two final-state momenta pi and p r with the help of an 
auxiliary light-like vector nf r [nf r = 0) using the usual Sudakov parametrization, 



Pi 



k 2 ^ 



P'r 



Z rPir + k±. 



k 2 ^ 



(3.1; 



where pf r is a light-like momentum that points towards the collinear direction and k± jT is 



the momentum component that is orthogonal to both pi r and rii r (pi r ■ k± t 



TU 



0). 



Momentum conservation requires that Z{ + z r — 1. The two-particle invariant mass of the 
collinear partons is 

k± r 

Sir = • (3-2) 



Zj Z r 



The collinear limit is defined by the uniform rescaling 

k i r ^ A/u i r i 



(3.3) 



and taking the limit A — > 0, when the one-loop squared matrix element of an (m + l)-parton 
process has the following asymptotic form [211]:^ 

2Re(M { Zl(PiiPn ■ ■ OI|A*£U(Pi,Pr, •••))- 



2e 



2Re(^)(p in ...)|P^ r |M«(p ir ,...)) 



where 



87ra s cr 



Sir J 



i r(i + e)r 2 (i- £ ) 



Pgj^(p ir ,---)) 



f(0) 



(3.4) 



(3.5) 



(4vr) 2 - £ T(l - 2e) 

The meaning of the ~ sign in Eq. (J3.4|) is that we have neglected subleading terms (in 
this case those that are less singular than 1/A 2 ). In order to simplify further discussion, 
following the notation of Ref. jSj, we introduce a symbolic operator C ir that performs the 
action of taking the collinear limit of the one- loop squared matrix element, keeping the 
leading singular term: 

C ir 2Re{M { V + i(Pi,Pr, ■ ■ .)\\MV +1 (jPi,Pr, ■■■)) = 

>\P}l\M%( Pir , 



87ra s /i 



2e 



$ i r 



•)) 



,2 \ e 



iira s c r ( — ) cos(Tre) (M < ^ ) (p tr , 



PZl\MM(p ir ,...)) 



. (3.6) 



' Since we deal with final-state singularities only, we have Si r > and we can write the usual factor 
(-/i 2 /s ir ) e as {pP/sirY cos(7re). 
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The m-parton matrix elements on the right-hand side of Eq. (|3.4J) are obtained from the 
(m + l)-parton matrix elements by removing partons i and r and replacing them with a 
single parton denoted by ir. The parton %r carries the quantum numbers of the pair % + r 
in the collinear limit: its momentum is pf r and its other quantum numbers (flavour, colour) 
are obtained according to the following rule: anything + gluon gives anything and quark 
+ antiquark gives gluon. The kernels Pf°J r and Pft r are the <i-dimensional Altarelli-Parisi 
splitting functions and their one-loop corrections, which depend on the momentum fractions 
of the decay products and on the relative transverse momentum of the pair. For the sake 
of simplicity, we label the momentum fractions belonging to a certain parton flavour with 
the corresponding label of the squared matrix element, Zf. = z*. In the case of splitting 
into a pair, only one momentum fraction is independent, yet, we find it more convenient 
to keep the functional dependence on both Zi and z r . Depending on the flavours of the 
splitting products the explicit functional forms are* 



{ f ,\P£l(z i ,z r ,k1;e)\u)=2C A 



1L 



klAPflX^ z r ,k1_;e)\v) = T R 



(r\P£Uz i ,z r ;e)\s)=8 rs C v 



-9 
l + zf 



Zi 



2(1 — e)zjZ r 



rv I rv I 



(3.7) 
(3.8) 
(3.9) 



where in the last equation we introduced our notation for the spin-averaged splitting func- 
tion, 



fifr 



\PfiU\ Z h z r, ^jj £ )) ■ 



(3.10) 



The one-loop kernels are 



(AilP&C*, zr, Aft e)\u) = 4% n ( Zi , z r ; e^P^Zi, z r , e)\v) 



„9i9r 

N.a 



1 2 ^ Z-) Zf 



h^ 1 h v 



-4C A r 

(A* I P^ilr ( Z i 1 Z r ) ^1 ' £ ) I U ) = r &%ni Z i 1 Z r ) I Pqiq r ( Z i 1 Z r 1 ^± i £ ) I V ) > 

(rlPjUte, z r] e)\s) = r** («, z r )(r\P£l(z t , z r ; e)\a) + 5 rs C F r% 



1 — ez r 



(3.11) 

(3.12) 
.(3.13) 



The r &iVon (zi, z r \ e) singular factors are expressed in terms of corresponding unrenormalized 
r & (zi, z r \ e) factors. The relation between the two forms is given by the equation: 



r a \ z i: Z r ]Sj 



A) s £ 



2e (47r) 2 c r L\ s 



2 \ £ 



COS! TIE 



(3.14) 



The unrenormalized r[^ r (zi, z r ) factors and the r{i( r factors can be trivially obtained from 
the T>j' ^f 'fr functions, that were computed in Ref . [2D] • in the case of gluon splitting those 



*We remind the reader that the formulae are valid in the CDR scheme. 
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functions are symmetric under the exchange of Zi and z r . To make this symmetry manifest, 
we have re-cast the original expression of Ref. [20] for the gluon splitting into a qq pair into 
an equivalent form which exhibits the z% <-> z r symmetry, 



1 



(Ca — 2C F 



m=l 



Li, 



Li, 



-3r 



2 £ 



11 4T R 2T R 

y C A — n f — n s - 3C F 



C A + 4T R (n f -n s ) 
3(3 - 2e) 



(3.15) 



In the case of gluon splitting into two gluons, the same symmetry is valid, which, however, 
we choose not to make manifest. Instead, we use a form, where the poly logarithms are 
regular in the z r — > limit (which will be convenient when we compute the soft limit of 
this expression in Sect. 13.1. 31 see Eq. ([3. 240 ). 



tie 



Sm(7T£ 



E 2£ 

m=l 



2m-l t ; 
ijl 



2m-\ 



(3.16) 



and similarly in the case of quark splitting, 



^ii ^r) 



tie 



sin ure 



+ J> m [(l + (-l) m )C, 



m=l 



2C F 



LL 



(3.17) 

Eq. (|3.16|) also shows that polylogarithms with even subscripts do not appear in the e 
expansion of the rl l9r singular function for gluon splitting. 

The r NS non-singular factors do not depend on the momentum fractions, 
C A (l-s)-2T K (n f -n s ) mr C A -C F 



„9i9r 



(l-2e)(2-2e)(3-2e) ' ' NS 1 - 2e ' (3 ' 18) 

In Eqs. (|3.15jl and (|3.18j) m and n s denote the number of fundamental fermions and scalars 
that can circulate in the loops. The case of QCD is obtained by setting n s = 0. 

The gluon-gluon and quark- ant iquark splittings are symmetric in the momentum frac- 
tions of the two decay products (even though this is not manifest for gluon-gluon splitting), 
while the quark-gluon splitting is not. Nevertheless, we do not distinguish the flavour ker- 
nels P qg and P gq . The ordering of the flavour indices and arguments of the Altarelli-Parisi 
kernels has no meaning in our notation, i.e., 



P fifr (zi, z r ] e) = P frfi (z r , z^ e) 



(3.19) 



Thus, it is sufficient to record the kernel belonging to one ordering. We keep this convention 
throughout. 
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3.1.2 Soft limit 



The soft limit is defined by parametrizing the soft momentum as p(f = Xqp and letting 
A — > at fixed qp. Neglecting terms less singular than 1/A 2 , it was found [2T] that 

2Re{M { V +1 (p r ,...)\\M% +1 ( Pr ,...))~ 



2= 



-87ra s cr 



i k=£i K 
1 



(0), 



..)\TiT h \MS(...)) 



C 



TVS 



A - n~. 7 s 

e z sm T£ 



COS 7T£ 



A) Se 

2e (47r) 2 c r 



i^i,k \ ' 



(0) / N|2 



l^k fc )(--0| 2 



(3.20) 



if r is a gluon. Similarly to the Cj r operator of taking the collinear limit, following Ref. jH] 
we introduce another symbolic operator S r that performs the action of taking the soft limit 
of the squared matrix element, keeping the leading singular terms. With this notation 
S r 2 Re(A4^ +1 (p ri . . OII-Mm+itjV, • • •)) 1S equal to the right hand side of Eq. ()3.20|) if r is a 
gluon and S r 2 Re(.M^+ 1 (p r , . . Oll^m+iCPn ...))= if r is a quark. 

In Eq. ()3.20)1 the m-parton matrix element on the right-hand side is obtained from the 
(m + l)-parton matrix element on the left-hand side by simply removing the soft parton. 
The eikonal factor is 

S lk (r) = —. (3.21) 

Si r S r k 

Note that Eq. (J3.20J) is valid only for the case of final-state partons. The general case can 
be found in 



3.1.3 Matching the collinear and soft limits 

If we want to regularize the squared matrix elements in all singly-unresolved regions of the 
phase space then we have to subtract all possible collinear and soft limits, i.e. subtract the 
sum 

E (E \ C * + S r ) 2Re(M^ + i(Pi,Pr, ■ ■ .)\\M% +1 ( Pi ,Pr, ■■■)), (3.22) 

where the 1/2 symmetry factor because in the summation each collinear configuration 
is taken into account twice. Subtracting Eq. ()3.22)) we perform a double subtraction in 
some regions of the phase space where the soft and collinear limits overlap. In order to 
compensate for the double subtraction, we need to find the collinear limit of the right 
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hand side of Eq. (|3.20J) when gluon r becomes simultaneously collinear to parton i. In 
deriving this limit we use that in the collinear limits (i) the factors multiplying two-parton 
colour-correlated squared matrix elements are independent of k, therefore using colour 
conservation (Eq. (|2.3jl ) we can perform the summation over k; (ii) the factors multiplying 
the three-parton colour-correlated squared matrix element are symmetric in k and / while 
l-M^L u n(f>i> • • -)| 2 * s antisymmetric, thus the sum of those terms is zero. Finally we have 



m;(i,k,l) ' 



C tr S r 2Re(M^ +1 ( Pi , Pr , . . .)\\M^ +1 ( Pi ,Pr, •••)) = 87ra s fi 2£ ——T. 

Sir %r 



2e 2 Zi ^ 



X 



2Re(^(^...)||A^W(^...)} 



-87ra s c r C A 



1 7T£ 



[J? Zi\ . . 8q S £ 
cos(iTe) + 



(•Maw--)! 1 



e*&nx\Tre) \s ir z r J ' ' 2e (47r) 2 cr 
Similarly, the soft limit of Eq. (|3.6|) when r is a gluon and z r — > is 



. (3.23) 



S r C ir 2Re(M^ + x(Pi,Pr, ■ ■ .)\\M% +1 (pi,p r , •••))= Sna s ^ £ — -T 2 



2 1 



x 



2Re(MS( Pi ,...)\\M$( Pi ,...)) 



-8na s c r C A 



1 



tie 



S, 



I COS{7T6) ~\ -r 

e 2 sin(7re) \s ir z r J 2e (Airj^cr 



|A<£W-0I 5 



. (3.24) 



Eqs. (|3.23|) and (|3.24|) differ by the term Zi = 1 — z r in the numerator of Eq. (J3.23|) . which 
is subleading if r is soft. Therefore, Eq. (|3.23|) can be used to account for the double 
subtraction: it cancels the soft subtraction in the collinear limit by construction, 



C ir (S r - C ir S r )2Re<^ 1 ||^ 1 ) = 
and the Cj r — Cj r S r difference is subleading in the soft limit, 

S r (C ir - C ir S r )2Re{M { Zi\\M% +1 ) = 



(3.25) 



(3.26) 



Accordingly, in order to remove the double subtraction from Eq. (|3.22|) . we have to add 
terms like that in Eq. (J3.23J) . That amounts to always take the collinear limit of the soft 
factorization formula rather than the reverse (like terms in Eq. (|3.24|) ). Thus the candidate 
for a subtraction term for regularizing the squared matrix element in all singly-unresolved 
limits is 



A 1 2Re(M { Zi\\M i ^ +1 ) 



(3.27) 



E? c ' 




2Re(M^ +1 (p i ,p r ,...)\\M^> +1 (p i ,p r ,...)) 



(i) 
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Note that the cancellation of the collinear terms in the soft limit actually requires the 
symmetry factor multiplying the collinear term, but not the collinear-soft one. This form of 
the Ai operator coincides with that derived in Ref. [S] for separating the singly-unresolved 
kinematical singularities of the squared matrix element at tree-level and is completely 
universal. 



3.2 Counterterms 



The expression given in Eq. (J3.27)) is defined only in the strict soft and/or collinear limits. 
In order to define true countertems, we have to extend it over the whole phase space. This 
extension requires an exact factorization of the m + 1 parton phase space into an m parton 
phase space times the phase space measure of the unresolved parton, 



d0m+i({p}) = d( fim({p}) [dp 



1 • 



(3.28) 



where we introduced the compact notations {p} = {pi, . . . ,p m +i} an d {p} = {pi, ■ ■ ■ ,pm}- 
Then the subtraction term that regularizes the kinematical singularities of the real- virtual 
cross section can symbolically be written as 



daZlf 1 = d0 m [d Pl ] Ai 2 Re(M 



(0) \\M {1) 

m+l\ K vl m+1 



(3.29) 



where we decompose the subtraction term as follows, 
A2Re(A«il|A<2i} = 



E 

r 

E 



E \ c t l 



(M)+^ 0,1) (M) 
E^ 0) (M)+f^)(M) 



E 



E C ^ 1,0) (M) 



(3.30) 



All terms above are functions of the original m + 1 momenta that enter the one-loop 
squared matrix element. The last terms in each line on the right hand side do not refer 
to the collinear limit of anything, but denote functions of the original momenta for which 
the notation inherits the operator structure of taking the various limits, but otherwise has 
nothing to do with taking limits. 

We now turn to the definition of each term in Eq. (J3.30j) . Each term will have the 
structure that a singular function (Altarelli-Parisi kernel or eikonal factor) is sandwiched 
between amplitudes. Both the singular functions and the squared matrix elements have 
their own loop expansions. The double superscipt on the subtraction terms refers to the 
number of loops in these loop expansions, the first one in the loop expansion of the singular 
factor while the second one in the expansion of the squared matrix element. 
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3.2.1 Collinear counterterms 



The collinear counterterms are 

Sir 



x2Re(A^W({pl^)|P)3 r (^ r ,^,^ !r ; £ )|^W({p}W)), (3.31) 
C£' 0) ({P}) = (Svra^^crcos^) 

S ir 

x(^H{p} (ir) )|PS r (^,^,^n £ )|A^S({p} W )). (3.32) 



The momentum fractions Zi r and 



and 2: rj j = — — , (3.33) 



V{ir)Q U{ir)Q 

while the transverse momentum k^ r is 

ijj, j- fi j- jU I j- ~fj, j- Vir j. Uir /r, q a 

_L ir >,i,rVr ^>r,iPi ' SirPi r j Si.r ^i,r 5 St - ,! *r,i • ( i 0.o4J 

CHi r y(ir)Q Ol ir y(i T )Q 

We used the abbreviations y ir = s ir /Q 2 = 2pi ■ p r /Q 2 , U(i r )Q = ViQ + VrQ with y iQ = 
2pi ■ Q/Q 2 , VrQ = %Pr ■ Q/Q 2 and is the total four- momentum of the incoming electron 
and positron, while pf r and a ir are defined below in Eqs. (J3.36J) and ()3.37|) respectively. This 
choice for the transverse momentum is exactly perpendicular to the parent momentum p^ r 
and ensures that in the collinear limit \\p%, the square of k^_ ir behaves as 

^Xir — $ir z r,iZi,r j (3.35) 

as required (independently of Q r ). In our computation the longitudinal component, propor- 
tional to Q r , does not contribute due to gauge invariance of the matrix elements, therefore, 
we may choose Q r = 0. The m momenta {p}^ r ' = {pi, ■ ■ ■ ,pi r , ■ ■ ■ ,p m +i} entering the 
matrix elements on the right hand side of Eqs. (j3.31|) and ()3.32|) are 

Pfr = rr— W+rf-airQ' 1 ), P% = T ^—P%, n^i,r, (3.36) 

where 

1 



V(ir)Q - ^y 2 ir )Q -^yir ■ (3-37) 



This momentum mapping leads to an exact factorization of the phase space in the form 
of Eq. (J3.28j) . The explicit expression for [dpi] reads 

[dPli(Pr,Pir] Q)] = jS(Pr,Pir; Q) ^ ^+(^r) » ( 3 ' 38 ) 
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where the Jacobian is 

Jl-m(Pr,PiriQ) = 



[1 - air) 



(m-l)(d-2)-l 



e(i - on 



2 \ l ~ yTrQ) a ir + Vrir + V*Q ~ Vr 

In this equation aj r has to be expressed as a function of the variable pi r 



(y r 



V 



irQ 



irQ 



Vr 



2(1 -VTrQ, 



(3.39) 



(3.40) 



It is straightforward to compute the e expansion of the collinear counterterms that we 
shall use later. We decompose the expansion into singular terms, finite contributions and 
terms that vanish as e — > 0, 

4°' 1} (W) = Voles C^\{p}) + KnCt l \{p}) + 0(e) . 



The pole parts can be written in the following unified form: 



Voles C£' 1] {{p}) = -87ra s /i 2£ — 

Si r 



(M^({p}^)\I({p}^;s) Pj^Ov, zr,i, fc Mr ; e) l-MgW* )) , (3.42) 



where§ 



J(M <,. ;£) ^ s ,(|)^( T ^ +7 ,I + g r , Ti I 



i/v 



with the usual flavour constants 



1 9 



00 



The poles of the C i7 ! ({p}) counterterm can be written as 



Voles C\l' Q \{p}) = -87ra sf i 2£ — ^S £ 



1 a s 



x 



(T 1 ^ _i_ / T 1 -2 /ti2 

I i r ii 

f / rri2 rri2 

-*■ i -*- r 



- i In y ir ^ + ^ ^ + 7r - 7tr) 



T?) In Zf. r +(T? 



T? + T? r ) In 2 r 



X {Ml\{p}^)\Pf fr { Zi ,ri *r,i, fc± Ar ; ^)I-M^({p} (lr) )) 



(3.41) 



(3.43) 



(3.44) 



(3.45) 

where T ir = T, + T r . After the cancellation of the poles is demonstrated (see Sect. 14.2)) . 
in a computer code one uses the finite parts of the counterterms. We collect all such finite 
contributions of the counterterms in Appendix [XJ 

§Note that I(m, e) = 1(e) + 0(e°) independently of m. 
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3.2.2 Soft-type counterterms 



We call the soft and soft-collinear counterterms soft-type terms because they all use the 
momentum mapping appropriate to the soft countertem. We define 

^(M) =-8 7 r« s/ i 2 ^^^(r)2Re(A^W({p}M)|T l T fc |^«({p}M)), (3.46) 



i k^i 



1 2z. 



C ir S^\{p}) = 8n as ^-^T^2Re(M^\{p}^\M^({p} ir) )) 



(3.47) 



i k^i 



1 TIE 

1 e 2 sinfTre) V 2 



Sik{r) cos(7T£:) + 



2s (47r) 2 c r 



1 718 ( 1 Zi r 



\M% >k) {{v} {r) )\ 2 



(3.48) 
(3.49) 



x 



C f 



Sm(7T£ \ Sir Z. 



cos(7re) + 



A) V- 2£ Se 

2e (47r) 2 c r 



The momentum fractions and eikonal functions were defined in Eqs. ()3.33j) and 1)3.21)1 . 
while the m momenta {p}^ = {pi, ■ ■ ■ ,p m +i} (Pr is absent) entering the matrix elements 
on the right hand sides of Eqs. ()3. 46)1 - 1)3.49)1 are defined as 



where 
and 



PC = A£[Q, (Q~Pr)/KM/K) , n ^ r , (3.50) 
K = a/1 -VrQ, (3-51) 



The matrix A£[K, K] generates a (proper) Lorentz transformation, provided K 2 = K 2 ^ 0. 
This momenum mapping leads to exact phase-space factorization in the form of Eq. ()3.28)1 , 
where 

[dpg»(Pr; Q)\ = rflfa Q) -j^hMpl) > ( 3 - 53 ) 
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with Jacobian 



(3.54) 



Finally, we record the pole parts of the soft-type counterterms: 
Voles Sj?V({p}) = Ma s ^Y,Y,\ S ^ r ) 



(Mi^{P}^\im (r) ;e),T l T k }\M^\{p}^)), (3.55) 



(3.56) 



2 Zj r ^2 



Voles C ir S^({p}) = -8na s ^--^Tl 

*{Mt\{P} {r) )\I{{P} {r) ;e) \Mt\{PY r) )), 
VolesS^\{p}) = 8™s/i 2e ^^k fc (r)|^^ fc) ({p}«)| 2 



2vr 



e 2 e y ir y kr e Q 2 J 2e 
2 zt 



Voles C ir S%°\{p}) = -Ma s ^-^T 2 \Mt\{p} {r) )\ 2 

Sir %r_i 



x^S £ 
2vr 



The finite parts are given in Appendix 1X1 



(3.57) 



(3.58) 



4 Counterterms for the integrated approximate cross 
section 



4.1 Factorization in the collinear and soft limits 

We wish to construct the approximate cross section ^ / x der^f 1 j Al by the same proce- 
dure we used to construct dcr^+'f" 1 , therefore, we start by studying the infrared limits of 

(M^ +1 \I(m, when the momenta of a pair of partons becomes collinear or when 

a gluon becomes soft. 
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4.1.1 Collinear limit 



In the limit when the momenta of partons i and r become collinear (as defined precisely in 
Eqs. dSU) and (Q) we find 1 

C ir {M^ +1 {pi,p r , ■ ■ -) \ I{m,e) \M^ +1 (pi,p r ,...)) = 87ra s /i 2e 
x- ((MSipir, ■ ■ -)l I(m,s) Pf} r \M<$( Pir , . . .)) 

+R ir (yir, Ziy {ir)Q , z r y (ir)Q ; m, e)(M^{p ir , . . .)| £j°/ r |M£?(Pin . . •))) ,(4-1) 

where the function R ir represents those terms that appear in addition to the usual collinear 
factorization formula of the squared matrix element, due to the presence of the insertion 
operator, 

«s ( At 2 



Ci(ziy {ir) Q;m,e) T- + C r (2r r y (ir .) Q ; m, e) T 2 - C (ir) (t/ ( , rW ;m,e) T 2 



+(T 2 ir - T- - T 2 ) S ir (y ir , ZiywQ, z r y( ir )Q; m, e) 



(4.2) 



The m parton matrix elements on the right hand side are obtained from the m + 1 parton 
matrix elements by removing partons i and r and replacing them with a single parton 
denoted by ir in the usual way. 

Note that the existence of a universal collinear factorization formula as given in Eq. (|4.1|) 
is not guaranteed by the factorization properties of QCD matrix elements, but depends 
also on the particular definition of the subtraction term da^'^ 1 , which determines the 
functional dependence of the insertion operator on the momenta. In being able to de- 
rive Eq. (J4.1j) it is crucially important that the part of I(m, e) that contains true colour- 
correlations, that is Su(yu,yiQ,yiQ;m,e), depends on its arguments only through the com- 
bination 

(4.3) 

yiQyiQ 

This expression is independent of the momentum fractions in the collinear limit £>i||p r - 
Consequently, the functions Su and S r / have the same limit as Pi and p r become collinear, 

CirSu = Cj r S r / . (4-4) 



^Note that the argument of the insertion operator on the right hand side of Eq. l|4.l[l is the same as the 
number of coloured external legs in the matrix clement. 
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This is important because coherent soft-gluon emission from unresolved partons implies 
that only the sum of \M^ liM \' + |^£- l!W) | 2 (or \M®. WA \* + \M^ r) \ a ) factorizes 



as 



(\M^ +im \ 2 + \M { Z W) \ 2 ) = f (MSlTtoTtPWjMS) . (4.5) 

Then, if and only if Eq. (J4.4p is fulfilled, we can combine the collinear limits as 

Cir (S,|-M^ 1;{ ,, | 2 + Sw|^^i ;W) | 2 ) = W^*,, j- {M^\T irTl Pf fr \M^) . 

(4.6) 

The insertion operators used in other general NLO schemes do not possess this property. 
4.1.2 Soft limit 

In computing the limit of (A4^ +1 \I(m, s)\A4^ +1 ) as the momentum of parton r becomes 
soft, we need the soft factorization formula for the colour-correlated tree amplitudes as can 
be found in Ref. [Bj. One finds 

S r (M ( ZM, ■ ■ 01 IK e) ...)) = -87roV e E \ S ^ r ) 



2 



><[(M^(...)\ 1 -{I(m,e),T t T k }\M^(...)} 



+^ikAyik^yir,ykr,yiQ,ykQ,yrQ;m J e)\M i ^ ii k) (. . .)| 2 ] (4.7) 

if r is a gluon and S r (A4^ +1 \I(m, s)\.Mm+i) = if r is a quark or antiquark. The m parton 
matrix elements on the right hand side are obtained from the original m + 1 parton matrix 
elements by simply removing parton r. In Eq. (J4.7|) the function 

ot s ( fi^ \ £ f 

RikAVik, Vir, Vkr, ViQ, VkQ, VrQ] m, s) = C A ^ S £ I —J [C g (y rQ ; m, e) 

+Sik(Vik, ViQ, VkQ] m, e) - S ir (y ir , y iQ , y rQ ; m, e) - S rk (y rkl y rQl y kQ ; m, e) ) (4.8) 



represents those terms that appear in addition to the usual soft factorization formula of 
the squared matrix element due to the presence of the insertion opeartor. 
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4.1.3 Matching the collinear and soft limits 



The collinear limit of the soft factorization formula Eq. (|4.7jl reads 

CfrS^gkfo.Pr, • • .)\I(m,e) \M { Zi(Pi,Pr, •••)) = 87ra s /i 2e — — T\ 

x ([Mt\Pi, ■ ■ 01 1(m,e) l-M^fe, •••)) + l-MffCft, . . .)\ 2 C ir R ik> ^ , (4.9) 

with 



CirRik,r(yik, Vir, Vkr, ViQ, VkQ, UrQ] m -, s) = ^A— S £ f — J 

x \p g {z r y{ir)Q] m, e) - S ir (y ir , Ziy {ir)Q , z r y {ir)Q ; m, e) ) . 
The soft limit of the collinear factorization formula Eq. (j4.1j) is 



(4.10) 



S r C ir (M { ^ +1 (pi,p r , . . .)| I(m, e) \M^ +1 (pi,p r , . . .)) = 8ira s fi 



(o) 



,2- 



2 1 



Sir %r 



x ((-M^fe, • • 01 1(jn,e) \M<$( Pi , ...)) + \M<$( Pi , . . OI 2 S r R*.) • (4.11) 



with 



S r R ir (y ir , Ziy( ir )Q, z r y^ r ) Q ; m, e) = S £ (j^J 

x (c g (z r y(i r ) Q ; m, e) - S ir (y ir , y(i r )Q, z r y(ir)Q] m, e 
Thus the same arguments as below Eq. ()3.24|) apply here as well, therefore, 

C ir (S r - Q3r)(ACil J K £ )l>fii> = > 
S r (C fr - C ir S r )(M^ + i\I(rn,e)\M^ + i) = 

and our candidate counterterm has the same structure as in Eq. (|3.27j) . 
Ai (M i £ ¥1 \I(m,e)\M^ 1 ) = 



(4.12) 



(4.13) 
(4.14) 



E 



J2\c ir + (s r -J2 c irSr) (\M^ +1 ( Pi ,p r ,...)\ 2 ^I(m,e)y (4.15) 



As before the cancellation of the collinear terms in the soft limit requires the symmetry 
factor multiplying the collinear term, but not the collinear-soft one. 
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4.2 Counterterms 



In order to define the counterterms, we extend Eq. (|4.15jl over the whole phase space as done 
in Sect. 13.21 We introduce the phase space factorization as in Eq. (|3.28j) and write the sub- 
traction term that regularizes the kinematical singularities of the integrated approximate 
cross section as 



(4.16) 



where 



■Ai(\m!% +1 \ 2 ® I{m,e 

r L i¥ zr 



(0,0®/) 
2 -ir 



(M) 
0) (M) 



«s(W(M)-E^^ w) (W) 



Sr m ({p})-J2C ir S^°\{ P })\ .(4.17) 

We now define all terms on the right hand side of this equation precisely. The structure of 
Eq. ()4.17|) is the same as that of Eq. ()3.30|) . Thus defining true subtraction terms starting 
from the limiting formulae of the previous subsection follows the steps of Sect. 13.21 



4.2.1 Collinear counterterms 



The collinear subtraction terms read 



4°'° 0/) (M) = 87ra syU 



2c 



x (M®m<V)\ I({p} [ir h m, e) P£i (z ijr , k ± ,y, e) \M<$({p}^)) (4.18) 



and 



C£ (M) = 8vra s /i 2e — Ri r (yir, z^ T y~ Q , z ryi y~ Q ; m, e) 

x(^£H^ (fr) )|^3 r (v.^.*Mr;e)l^£H{^ (ir) )>- (4-i9) 

Computing the pole parts of Eqs. ()4.18jl and ()4.19j) . we can easily see that 



and 



Voles 



Voles 



c£ ml \{p}) + c^\{p}) 
Rxm ({ P }) + c^\{p}) 



c 







(4.20) 



(4.21) 



Therefore, the sum of Eqs. (pHT^) and (f4~T8j) as well as that of Eqs. (f3~32jl and <$JM is finite 
in d = 4 dimensions. 
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4.2.2 Soft-type counterterms 



The soft-type counterterms are defined as 

s<^(M) = -W ta ££W r ) 



i k^i 



X 



(^(^^{^{P}^^^),^}!^)^)), (4.22) 



1 2* 



2 £ - 1 - ^i,r 



Sir %r,i 

xT? (-M^({p}M)|I({p}«;m,e)|>lW({p}M)), 



X Rik,r(Vik, Vir, Vkr, ViQ, VhQ, VrQ] m, e) 



h\ 2 



5 7.7* 



2 \ e 



X [Cg(z rt iy^ r )Q;m,e) — Sir(yir, Zi,rV(ir)Q, Zr,iy(ir)Q',™, 



(4.23) 



(4.24) 



(4.25) 



We can simplify the arguments in the second line of Eq. (|4.25J) using Eq. ()3.33|) . 

C g (z r>i y(i r )Q; m, e) - S ir (y ir , z i:r y( ir ) Q , z^y^Q; m, e) 
= CgiVrQ] m, e) - S ir (y ir , y iQ , y rQ \ m, e) . 

Similarly to the collinear cases, the pole parts cancel term by term between Eqs 
(l3~47H) and Eqs. KTR - KZty . 



Voles 



Voles 



Voles 



Voles 



o, 



^x(o,o) (M)+(S (i,o )(M) 



C ir S^({p})+C ir S^({ P }) 



0. 



(4.26) 

(4.27) 
(4.28) 
(4.29) 
(4.30) 



Consequently, the sum of Eqs. (g~g)}) and lf£22|) . that of Eqs. (l3~47j) and fl03J), that of 
Eqs. fEHH]) and f|4~^j) and that of Eqs. (J3~4T?|) and are finite in d = 4 dimensions. 

We help the reader in grasping the various cancellations by visualizing the subtraction 
terms in Fig.^ The first picture corresponds to the squared matrix element of the m + 
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Figure 1: Graphical representation of the squared matrix element and the subtraction 
terms. 



1 final-state partons. The fully shaded circle represents the Born amplitude, while the 
circle with a whole is the one-loop amplitude. The following terms in the first squared 
brackets correspond to the terms that build A-i 2 Re(Ai ( ^ ) +1 \\Ai^ +1 } J defined in Eq. fl3.30|) ; 

the first one represents (0,l)-type terms like C^'^iip}) and the second is for (l,0)-type terms 
such as Cfr ({p}). The first picture in the second line is the result of the integration in 
Eq. (|2.21|) . Finally, the terms in the second squared brackets represent terms that contribute 

to A-i (\-Mm+i\ 2 ® I{ra, e) J , defined in Eq. (14. 17|) : the first representing terms of the (0,0® 

J)-type, such as C^'° ({p}), while the second one stands for R x (0,0)-type terms like 
C 4 y X ^ '°^({p}). The factorized one-parton factors correspond to Altarelli-Parisi kernels, with 
azimuthal correlations included, or eikonal factors, with colour-correlations included, either 
at tree-level (fully shaded circles), or at one-loop (circles with holes), or i?-functions (boxes). 



The cancellations of the e poles occurs vertically, term by term as shown previously. 
The regularization of the kinematical singularities takes place horizontally, between the 
first term and the terms in the following brackets, separately in each line. Kinematical 
singularities, introduced by the subtractions terms, are screened by the jet function, J m , 
just as in NLO computations. 

We conclude that the (m + l)-parton contribution in Eq. ()1.4|) is free of e poles as well 
as unscreened kinematical singularities. We can set e = and compute the integral of 
do"^ 1 ^ by standard Monte Carlo methods. Such an integration uses the finite parts in 
the e-expansion of the subtraction terms as e — > 0. The finite parts of the integrated 
approximate cross section of Eq. (j2.21|) can be found in Ref. [H]," while those of Eq. ()3.30|) 
together with the finite parts of the counterterms in Eq. f)4.17|) are presented in Appendix 1X1 

"Notice the shift of m by one unit in this paper (m — > m + 1) as compared to the value in Ref. 
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5 Numerical results 



In this paper we have defined explicitly all subtraction terms that are needed to make da^_ 1 
integrable in d = 4 dimensions for processes without coloured partons in the initial state. 
We have proven the cancellation of the e poles explicitly. In order to demonstrate that 
the subtraction terms indeed regularize the cross section for the real-virtual correction, 
we consider the contribution to the theoretical predictions for the three-jet event-shape 
distributions thrust (T) and C-parameter in electron-positron annihilation, when the jet 
function is a functional 

J n (pi, ...,p n ;0) = 5(0 - 3 (pi, . . . ,p n )) , (5.1) 

with 3 (pi, . . . ,p n ) being the value of either r = 1 — T or C for a given event (pi, . . . ,p n )- 
Starting from randomly chosen phase space points and approaching the singly-unresolved 
(soft and/or collinear) regions of the phase space in successive steps, we have checked 
numerically that the sum of the subtraction terms has the same limits (up to integrable 
square-root singularities) as the squared matrix element itself. 

The perturbative expansion of the n th moment of a three-jet observable at a fixed scale 
ji = Q and NNLO accuracy can be parametrised as 



<0»)= [dOO»~(0) 
J cr Q dO 



^) A® + (^) 2 fl? + (^) 3 <#> , (5.2) 
where according to Eq. (|1.2jl . the NNLO correction is a sum of three contributions, 

Mn) _ n {n) n (n) n {n) ,~ „\ 

Carrying out the phase space integrations in Eq. ()1.4|) . we computed the four-parton con- 
tribution Cq)^ as defined in this article. The predictions for the first three moments of r 
and the C-parameter, obtained using about four million Monte Carlo events, are presented 
in Table ^ In performing the numerical integrations, we do not encounter more severe 
numerical problems than known from computing the real-emission contribution in NLO 
computations and the computation of differential distributions does not pose any problem. 
The required CPU time is however much longer because of the much more cumbersome 
expressions that enter the various loop matrix elements. 



6 Conclusions 



In a companion paper ^H] we set up a subtraction scheme for computing NNLO corrections 
to QCD jet cross sections to processes without coloured partons in the initial state. The 
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Table 1: The moments C^l and C ( 



n 


W;4 


r<( n ) 

°C;4 


1 


-(1.23 ±0.01) 


10 3 


-(4.33 ±0.05) 


10 3 


2 


-(2.55 ±0.01) 


10 2 


-(3.25 ±0.02) 


10 3 


3 


-(4.79 ±0.03) 


10 1 


-(1.80 ±0.01) 


10 3 



scheme is completely general in the sense that any number of massless coloured final-state 
partons (massive vector bosons are assumed to decay into massless fermions) are allowed 
provided the necessary squared matrix elements are known. 

Three types of corrections contribute to the NNLO corrections: the doubly-real, the 
real-virtual and the doubly-virtual ones. In Ref. jTUj we defined the subtraction terms for 
the doubly-real emissions; those to the real-virtual correction can be found in the present 
paper. By rendering these two contributions finite in d = 4 dimensions, the KLN theorem 
ensures that for infrared safe observables adding the subtractions above to the doubly- 
virtual correction that becomes also finite in d = 4 dimensions. 

The subtraction terms for the real-virtual corrections presented here are local in d = 
A—2e dimensions, include complete colour and azimuthal correlations. The expressions were 
derived by extending the singly-unresolved limits of the one-loop squared matrix elements 
over the whole phase space and also extending the singly-unresolved limits of the integrated 
approximate cross section, used for regularizing the kinematical singularities of the cross 
section for doubly-real emmissions over the singly- unresolved regions of the phase space. 

In order to demonstrate that the subtracted cross section is indeed integrable, we have 
computed the corresponding contributions to the first three moments of two three-jet 
event-shape observables, the thrust and the C-parameter. In performing the numerical 
integrations, we do not encounter more severe numerical problems than known from NLO 
computations. The required CPU time is however much longer because of the much more 
cumbersome expressions that enter the various loop matrix elements. 
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A Finite parts of the subtraction terms 



In this appendix we present the finite parts, as defined in the decomposition (j3.41j) . of 
Eqs. f!3.30|) and f)4.17|) . term by term. We start with the collinear counterterms. The finite 
part of the sum of Eqs. (j3.31|) and (j4.18|) is 

. o 1 



Tn 



C i; (M) + c™(M) =4 

2 MMSm^lPffir&r, Zr t u Ki, r ) Tin \M^{{p} {ir) )) 

+ Yl ( nn Sik (y~i~ki y~iQ> vug* m ) + cs ( m ) 



x {Mt\{PY lT) )\T,T k Pfl{z^ z r „ k ±>i>r ) \M^({p}^)) 



-Tin 



Ci(y~iQ; m ) + CS(m) 



X T\ (M%mW)\ Pftfar, Zr,i, fex Ar ) \M^({p}^)) 



(a.i; 



where the finite part of the one-loop amplitude is defined by the following equation: 

\M$({p})) = -U({phe)\M%({p})) + ^Tn \Mg\{p})} , (A.2) 

with I({p}] e) given in Eq. ()3.43|) . The finite parts of the functions + CS and C, + CS 
can be found in the appendix of Ref. [TT]. Next we turn to the finite parts of Eqs. ()3.32|) 
and (J4.19|) . which are given separately for the various flavour combinations. For fi = f r = g: 

1 



Tin 



c^({p}) + & 0) ({p}) 



9i9r 



C A Aott 



x 



(^Fin X - (c gg {z^y TrQ - m) + C gg {z r>i y^ Q ; m) - C gg (y~ Q ; m) 
+ n f ( G qq( z i,rV^ Q ; m ) + C ! g? («r,i2/fr q) m) - C q g{y^ Q ; m)J 

— ^irKViri z i,rV'^ r Q^ ^t0^q\ m ) ~ CS(m) 



1 , D Zi 



--ln 2 ^ 

2 z r i 



In 2 y ir - \Tvy ir ln(z i>r z r> i) 



vr 2 I 
3 + 2C A Q 2 



-|(C A - 2T R n f ) (M%mW)\ W({p}W)> 



(A.3) 
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For fi = g, / r = g: 



cS£ } (M) + « u,u; ({p}) 



Kx(0,0), 



4cC — 



x [ C F Tin C qg (z i>r y^: Q ; m) + C qg (z r ^ Q ; m) 

-C A .Hn ^C 99 (y~ Q ; m) + n f C 9g -(^ Q ; m) 
+(C A - 2C F ) .Hn Sfr^i,., Zi, r y~ Q , Zr,iy~ Q ; m) + CS(m) 
+(C A -2C F ) Qln 2 y ir -y + 1^ +2( lg - lq )(2-\ny ir ) 

(\ Z ir 7T 2 \ 1 

-C A ( - In 2 j£ + lny ir ln(z ijr + — J + -(C A + 4T R n f ) 



(A.4) 



For fi = q, f r = g: 



Tin 



c^({ P })+cZ m ({p}) 



4a 2 — 

Si r 



Cf Tin 

-C\ Tn 
-C\ Tn 



C qg {z i>r y^ Q ; m) - C qg (y TrQ ; m) 



+2C F Li 2 



$>i r (yi r , Zi jr y ir Q, z r ^y ir g] m) + CS(m) 
- In y ir In z % 



Zr.i 



Zi,., 

( \ In 2 ^ + \ In 2 y ir + In ^ In + 2 Li 2 f - 

+yln^)(^)({p}^)|^)(^,, rii ) |M£)({p}^)> 



7T _ 

~3 



+c F (c A -c F )i^)({^))r 



(A.5) 



The case of fi = g, f r = q is obtained from the latter with the i <-> r substitution. 
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The finite parts of the soft subtractions are as follows: 



Tin 



x 2He<^g)({p}W)|T i T fc J?n|Ml ) ({p} (r) )> 



+ E ( [ S iK^^ Q ,%;^) + CS(m)]|^ fe)W) ({p}W)| 2 



+Tin 



C i (y jQ ;m) + CS(m) T? |^ (iiJfc) ({p}<*' a 



(A.6) 



Tin 



4 1,0) (M) + «C ( °" 0) (M)J = 4a s 2 E E ^ ^(r) |^S i>fe) ({p} ( ^' 2 



, 1 , 2 Vik 7T A) i A* 
x | - In — m — 

2 y ir y fcr 3 2C A 



Tin 



; C gg(yrQ; m) + n{C q q(y rQ ] m) 



-Tin \^>ik{yik, UiQ, VkQ\ m) - S ir (y ir , y iQ , y rQ ; m) - S kr (y kr , y kQ , y rQ ; m) - CS(m) 



-E^i^S,,-,(w M )i 2 } 



(A.7) 



Finally, we present the finite parts of the collinear-soft subtractions: 
Tin 



c ir sr\{p})+c ir sr^\{ P }) 

2Re{M%\{p}^)\Tin\M£\{P}^)) 



Sir %r 



X 



+ E ( E nn [Ml®, l/fo, ifo,; m) + CS(m)J |A^ - ({p} (r) ^ 2 



+Tin 



C^m) + CS(m)l T* |A<£>({p} (r) )l 2 , 



(A.8) 
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Tin 



C ir S^\{ P })+C ir S^\{p})\ =-W s C A -^T*\Mt\{PY r) )\ 2 

Sir 



x - In 



1 , 2 ^i,r 



m — - — Tin 



3 2C A Q 2 



-C 99 (2/ rQ ; m) + n f C q q(y rQ ; 



m 



+Tn 



SiriVir, ViQ, DrQ] ™) + CS(m) 



(A.9) 



For the sake of completeness, we recall the finite part of the first two terms in Eq. (J 1.4)) 
from Ref. adapted to the present case: 



L uo m+l + ucr m+l 



^Ke)] =-A/"ip V d0 m+ i— - — 



{m+1} 



x<!2Re(^ + i(W)l^l-M m+ i(W)) 



(i) 



^Tn [S ik (y ik ,y iQ ,y kQ ;m) + CS(m)] |A^2 + i ;(iifc) (M)| 



-Jin [C J (y iQ ;m) + CS(m)] T 2 \M { °U{p})\ 



(A.IO) 
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